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Abstract 

A new gauge fixing condition is discussed, which is (lattice) rotation invariant, has the 
'smoothness' properties of the Landau gauge but can be efficiently computed and is 
unambiguous for almost all lattice gauge field configurations. 
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Introduction 

In gauge theories gauge fixing is necessary to make the perturbation expansion well 
defined. One of the advantages of the lattice regularization of gauge theories, compared 
to the continuum formulation, is that gauge fixing on the lattice is not necessary. This 
is true if only gauge invariant operators are considered. In a growing number of cases, 
however, it is found advantageous to also use gauge variant operators which makes it 
necessary to fix the gauge even on the lattice. 

For instance, it is interesting to study the gauge variant propagators of quarks and 
gluons [|^]. Also when one wants to compare with gauge dependent results from pertur- 
bation theory gauge fixing on the lattice is necessary. In other applications one leaves 
out the strings of gauge parallel transporters that make operators gauge invariant, like 
in point split fermion bilinears or extended glueball operators. Since this eliminates the 
fiuctuations of the gauge fields contained in the string, it significantly improves the signal 
to noise ratio in a Monte Carlo experiment 0. To restore gauge invariance one must 
now fix to some uniquely defined gauge. Popular choices are the Landau or Coulomb 
gauge. An alternative picture justifying this strategy of leaving out parallel transporters 
is that by fixing to the Coulomb gauge, one dresses the bare charge with its physical 
gauge boson cloud instead of squeezing this cloud through a connecting fiux string. For 
noncompact QED this interpretation of the Coulomb gauge in terms of a photon cloud 
can be substantiated ^, ^ . 

Let us consider the non-gauge fixed gauge field U^{x) G G with gauge group G = 
SU{n) or U{n). Suppose the gauge fixed configuration minimizes a (gauge dependent) 
function R{U). An example is the lattice Landau gauge fixing condition, where 

R{U) = ReJ2[l--tTiU,ix))]. (1) 

X,IJ. '^ 

The gauge fixing algorithm should provide the transformation Q such that R{U^) as- 
sumes its minimal value, where 

U^ix) = n{x)U^ixMx + iiy. (2) 

If the starting configuration is gauge transformed with g, it should find the transformed 

n%x) = gon{x)g{xy, Q^^xY = g{x)n^x)gl (3) 

where the global gauge transformation gQ need not be fixed. The algorithm should find a 
unique solution. This means that for all gauge equivalent initial configurations, it should 
find the same gauge fixed configuration (up to a global gauge transformation f^o)- It is 
clear that this is impossible in principle if the absolute minimum of R is degenerate. 
But if this would be the only source of ambiguity, one would argue that such an exact 
degeneracy is an exceptional situation which will not occur in practice. However, finding 
the absolute minimum of R is extremely difficult and iterative minimalization procedures 
usually run into some local minimum. Finding the absolute minimum is a numerical 
problem of non-polynomial complexity. 
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In the continuum limit, the local minima of R correspond to solutions of the dif- 
ferential Landau condition and the existence of multiple solutions was first discussed 
by Gribov |^. Also on the lattice we shall call gauge fields for which R has a (local) 
minimum Gribov copies. The existence of Gribov copies on the lattice is demonstrated 

In order to remove the Gribov ambiguity it has been proposed to first fix to a uniquely 
defined axial gauge Pj. But even if the starting gauge field is unique, one can still run 
into different local minima by using different gauge fixing algorithms. Moreover the axial 
gauge breaks rotation invariance, and some effect of this might be present in the local 
minimum that is favored by this particular starting configuration. 

Another approach to remove the ambiguity is to include in the gauge fixing condition 
the requirement of averaging over all initial gauge transformations, such that one aver- 
ages over a large number of local minima. In practice this procedure is too expensive. 
Furthermore it is to be expected that the set of local minima that will be found depends 
on the algorithm used, e.g. it may be different for local algorithms and for non-local 
ones like multigrid 0. Finally, averaging over many local minima re- introduces noise in 
gauge variant operators 0. 

In the following we shall discuss a gauge fixing condition which can be computed 
efficiently and unambiguously for all but a set of exceptional configurations. It shares 
desirabel properties, like smoothness of the gauge fixed configuration, with the Landau- 
like gauges. Of course, it should be noted that we are not solving the Gribov problem 
for the usual, local. Landau and Coulomb gauge. We avoid the problem by constructing 
a modified gauge condition which is no longer local. Because of this non-locality it is 
not clear if our gauge fixing condition is well suited for perturbative calculations. If not, 
this will complicate a comparison of gauge dependent results from perturbation theory 
with nonperturbative lattice calculations. Unlike axial gauges our gauge fixing does not 
violate global symmetries like e.g. rotation invariance. Therefore it is well suited for 
the construction of improved interpolating fields for hadrons. To project on a definite 
symmetry sector of the physical Hilbert space it is essential that the hadron operators 
transform appropriately under the various global symmetries. It should be noted that 
the transfer matrix formalism requires the operators to be well localized in euclidean 
time. This suggests to fix to Coulomb- (rather than Landau-) like gauges when one has 
this application in mind. 

Laplacian gauge fixing 

For a given gauge field configuration, we want to find a gauge transformation Q E G 
which transforms as in (^ and captures the physical 'smoothness' of the gauge config- 
urations. Relaxing the property that Q E G, the above requirements are reminiscent 
of finding covariant renormalization group transformations or multigrid blocking kernels 
and this may suggest to consider eigenfunctions of the gauge covariant laplacian A{U) 
as in []TU|, [TT], |T2|. The covariant laplacian (up to a minus sign) is defined by 

A{U){x, yr" := J2mx - y)6''' - U,{xr'6{x + fi - y) ~ U.iyy'^'Siy + fi - x)] (4) 



and it is a positive (semi) definite operator. We consider a finite lattice with lattice 
spacing a = 1 and with periodic boundary conditions. The eigenfunctions /'' of A(^) 
are defined by 

AiU)ix,yrfiy)l = yfix)l (5) 

where the A* > are the eigenvalues. We shall often suppress a summation over repeated 
indices. For a non-degenerate spectrum, these eigenfunctions are unique up to a global 
phase, 

fix): ^ eM^vlfix):. (6) 

Under gauge transformations U —>■ U^ the eigenfunctions transform as f{x): — *■ 
g{xY^f{x)l. This implies that the GL{n) matrix field 

M{xr := f{x)t (7) 

has the same transformation behavior as gauge transformations Q{xy°''', cf.(^. The 
index b in (|^) labels a set of n eigenfunctions. The transformation behavior of M does 
not depend on which n eigenfunctions we choose, but in order to select the smooth modes 
in the gauge field, we shall choose the eigenfunctions with the lowest n eigenvalues. This 
is in the spirit of refs. |jlO|, |ll], |^ and will be discussed further below. 

In general the matrices M(x)"'' are not yet in the gauge group G. However, almost 
all M can be uniquely decomposed as M = Q^R with fi G G and transforming as in (^. 
To see this for the case of G = SU{n), we first write M = WP with W a unitary matrix 
and P a positive hermitian matrix. This decomposition is unique if 

P = (M^M)'/' = 5tdiag(/if , ...,fil/')S (8) 

is invertible, i.e. if M has a non-zero determinant. The unitary matrix S in (H) diago- 
nalizes M'^M and the Ha are the eigenvalues. Then W = MP~^ and Q^ = W exp{—ia) 
with expm = diag(det(iy)^/"', ..., det(iy)^/"). If the gauge group is U{n) the determi- 
nant need not be factored out. 

More explicitly we define the gauge transformation Q leading to our Laplacian gauge 
by 

Q^{xy^'' = f{x)l{P{x)-y\-'''^''\ (9) 

It is clear from the construction that this fi^ indeed transforms as in (^), whereas P and 
a are invariant. 

The global U{1) transformations of the eigenfunctions (|^) induce global abelian trans- 
formations of f^A. This follows from the observation that M"^ -^ M'*^exp(i?7^) implies 
that S""^ — >■ exp{—i'r]°')S"''' exp{ir]^) and consequently P"^ -^ exp{—iri"')P"'^ exp^irj'') as 
well. From (j^) it then follows that 

n^{xy' -^ exp(-zr/" - fi)nAixr\ (10) 

where rj = J2a v"'/n. 

In general there is no global SU{n) invariance, only if the lowest n eigenvalues are 
degenerate one can show along similar lines as above, that a global SU{n) invariance is 
present. 



Singularities and ambiguities 

The construction of Q/^ outlined above is unambigous for almost all gauge field con- 
figurations: For a given gauge field U the spectrum of A{U) is invariant under gauge 
transformations g and if the spectrum is non-degenerate, the eigenvalues provide an 
unambiguous labeling of the eigenfunctions from which the ^a's are constructed. In 
particular, the n eigenfunctions with the lowest eigenvalues can then be computed unam- 
biguously (up to a global phase). The subsequent decomposition works if det(M(x)) 7^ 
for all X. 

It is clear form this resume that there are two sources of ambiguity in our construction. 
If the (ra+l)-th eigenvalue is degenerate with the n-th (labeling them in ascending order), 
the prescription to take the 'lowest ra' is not well-defined. This means that Q/^{U) 
becomes ambiguous at those points in configuration space, where level n and n + 1 cross. 
The level crossing condition fixes one degree of freedom. Hence the singular manifold 
is of co-dimension one. The second source of ambiguity is at points in configuration 
space where one or more of the det(M (t/)) become zero. The condition that the complex 
determinant vanishes at some point x fixes two degrees of freedom and hence the singular 
manifold is of co-dimension two. In addition, for SU{n) the determinant has to be 
factored out. This requires to take the n-th root of it, which is well-defined only if the 
determinant is not negative real. This excludes a singular manifold of co-dimension one. 
In the continuum limit it is expected that these singular sub-manifolds become Gribov 
horizons giving configuration space the required non-trivial topology to incorporate flux 
sectors and 6 vacua [r3|. 



This discussion shows that the probability to find singular configurations in Monte 
Carlo simulations on a finite lattice is zero and our gauge fixing condition is unambiguous 
for all practical purposes. 

Relation to the Landau gauge 

It is an important property of the Landau gauge that the gauge fixed configuration is 
smooth in the sense that the link field U^j,{x) has similar fluctuations as the gauge invari- 
ant plaquette field. Of course there will be localized fluctuations as well as obstructions 
due to nontrivial topology (instantons), but on the average U is close to one if the average 
plaquette is close to one. By constructing gauge transformations from the eigenfunctions 
of the covariant laplacian and selecting the modes with smallest eigenvalues, one expects 
that the gauge fixed configurations are again smooth in the above sense. 

To illustrate this we shall first consider pure gauge configurations, U^{x) = g{x)g{x + 
fiy . Here the Landau gauge fixing condition has the solution U^{x) = 1. The eigen- 
functions of the laplacian are gauge transforms of plane waves, 

fix): = g{xr'exp{tkx), 

X^ = 5:(2sin(V2))^ (11) 

where the lattice volume is V^ = A^'^ and the momenta are k^ = 2'Kn^/N , n^ = 0, ..., N—1. 
The label s is such that s = 1, ...,n corresponds to the n zero momentum modes. For 



the pure gauge case the construction of Qa is simple and one finds 

Q^{x) = QogixY, (12) 

where Qq is a global SU{n) transformation which is not specified by our construction. 
Note that these are precisely the gauge transformations that fix U to the Landau gauge. 

Since the space dependence of the eigenfunctions for this simple case comes in as a 
t/(l) factor which drops out in the reduction to SU{n) we cannot illustrate the impor- 
tance of taking the low eigenvalue modes here. With g G U{n), however, our construction 
gives fl^ix) = Qog{xy exp{ikx), if eigenfunctions with momentum label k are used. Only 
ioT k = this corresponds to the Landau gauge, for |A;| > the gauge fixed configuration 
is less smooth. 

To elucidate the generic case, we consider the function 

x,y 

Minimalizing E under the constraint that "^Zx f {^Ya f i^^a ~ V^'^* provides the eigen- 
functions of A(f/) as solutions (normalized to the volume V). After substituting the 
decomposition (|^) in (|TB[) one finds 

y\^ = 2nRey [1 - -ii{e'''^''^ P{x)n^{x)UJx)^r^{x + ^)^P(a; + /i)e-^°(^+'^))] (14) 

— ^ — ^ Tt 

a x,fi '" 

This leads one to consider a modified Landau function i?A, 

i?A(f/") =ReY,[l--ii{W^{x)n{x)U^{x)n{x + fi)% (15) 

The weight function W^{x) is gauge invariant and follows from the eigenfunction con- 
struction, 

W^{x) = P{x + /i)P(x)e'("(")-"("+^)\ (16) 

The minimum of R^ over gauge orbits fi is assumed for Q = Q^ , coming from the eigen- 
function construction: i?^*" = J2a=i ^af^n. The weight function Wfj,{x) is a function 
of all gauge links Ufj,{x) and in this sense our gauge condition is non-local. Clearly R^ 
reduces to the (local) Landau function (0) if Wfj,{x)ab — > Sab- 

This shows that our laplacian gauge condition corresponds to finding the absolute 
minimum of a modified Landau function (p!5|). Trying to find the absolute minimum of 
i?A for a fixed W^{x) is practically impossible, but by using the W^{x) which follows 
from the eigenfunction construction, the problem is easily solved. 

Discussion 

We have shown that a smooth gauge, similar to the Landau gauge, can be defined 
unambiguously on the lattice, except for configurations which have measure zero in 
configuration space. In order to find it for gauge group SU{n) or U{n), one needs to 



compute the eigenf unctions corresponding to the n smallest eigenvalues of the covariant 
lattice laplacian and make the reduction (^) to SU{n) or U{n) matrices. The time 
consuming part is the computation of the n lowest eigenf unctions. They can be computed 
with the Lanczos algorithm, or they can be found by inverse iteration (see e.g [|12]), after 
the n eigenvalues have been computed with the Lanczos algorithm. Since the laplacian 
is hermitian, this should be feasible even on large lattices. 

The main virtue of our gauge fixing procedure is the possibility to uniquely specify 
a smooth Landau-like gauge in non-perturbative lattice computations without breaking 
global symmetries like rotation invariance. At present it is not clear (to us) if laplacian 
gauge fixing can also be used in perturbation theory. Even if this is not possible, laplacian 
gauge fixed fields may be useful in the construction of improved hadron or glueball 
operators. For this purpose the Wuppertal group 0] has used the inverse covariant 
laplacian before, which is reminiscent of our approach. A construction based on laplacian 
gauge fixing, however, offers the possibility to interpret results physically in terms of 
gauge fixed quark and gluon fields. 
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